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ANALYSIS OF THREE LIQUID-METAL QUALITY MEASUREMENT TECHNIQUES 

by Ted W. Nyland 
Lewis Research Center 

SUMMARY 

Three techniques were theoretically evaluated to determine potential accuracy and 
design problem areas. The techniques considered were calorimetric, flow measure- 
ment, and acoustic techniques. Calorimetric techniques had root-sum-square errors 
ranging from 0. 0028 to 0. 015 for heating and condensing calorimeters, to 0. 0027 for 
throttling calorimeters. These estimates were based on applying systematic corrections 
for heat losses and did not include sampling errors. Sampling errors may approach 0. 01 
or more. Two-phase flow measurement techniques have potential root-sum-square er- 
rors of 0. 014 in a small range of weight flows around the design point. Problems asso- 
ciated with achieving thermodynamic equilibrium during the passage of an acoustic pulse 
in a two-phase fluid apparently limit the application of sonic velocity measurement tech- 
niques. 


INTRODUCTION 

Measuring the relative composition of two-phase fluids is a problem that has chal- 
lenged instrument designers for years. This problem has been encountered in a wide 
variety of technologies and has resulted in measuring techniques suitable for environ- 
ments from steam power systems to cryogenic flow lines (refs. 1 to 9). One of the most 
difficult environments yet encountered is that of Rankine-cycle space power systems. In 
this case the fluids are alkali metals, most commonly potassium and sodium, and the 
temperatures range from 1500° to 2200° F. The literature cites few cases in which the 
quality or liquid-vapor content of such fluid streams is measured. An exception is ref- 
erence 10, in which measurements with a throttling calorimeter are reported. How- 
ever, details on the instrument design and estimated accuracy are not included. Instru- 
ment design and error analyses of quality measurement techniques are reported in ref- 
erences 11 and 12. More commonly, liquid- vapor content is estimated from overall 


heat-balance calculations, a procedure limited in accuracy by extraneous heat losses. 

Quality measurement is vital to the Rankine-cycle space-power-system program be- 
cause it defines the energy content of the two-phase fluid. Space-power-system studies 
(ref. 13) indicate that operation with turbine-inlet vapor of very high quality to a slightly 
superheated state is required to achieve a minimum system specific weight. For system 
testing with wet vapor, quality measurements are required to determine the turbine-inlet 
conditions and to monitor boiler performance. Quality measurement is also an aid in 
boiler development work because it provides a measure of boiler performance. 

A quality meter is not sufficient to fill all the needs for two-phase instrumentation. 
Instrumentation to define flow regimes (i. e. , bubbly, slug, annular, or mist flow) and 
void fraction (percentage vapor by volume) is required in studies of the boiling and con- 
densing processes. Also, instrumentation to detect liquid droplets in vapor streams and 
to measure droplet sizes and velocities is of interest, especially in turbine -blade - 
erosion studies. Quality measurement does, however, provide a most basic need in two- 
phase instrumentation, which is energy-content measurement. 

This report analyzes three measuring techniques that are applicable to the measure- 
ment of quality in alkali-metal two-phase streams, (1) calorimetric, (2) weight-flow 
measurement, and (3) sonic velocity or acoustic techniques. Quality is the ratio of the 
weight flow of vapor to the total weight flow of the stream. One prime prerequisite for a 
liquid- metal-quality measuring system is that it be inherently capable of measurement 
without the need for calibration. Each technique is analyzed from the standpoint of the 
fundamental principles involved, the physical measurements required, and the practical 
problems of applying the technique to two-phase alkali-metal systems. 

For the calorimetric and two-phase weight-flow techniques, error equations are de- 
veloped, and estimated probable errors are calculated. For the sonic velocity tech- 
nique, the basic assumptions required for quality determination are examined in relation 
to the properties of the fluid stream. In this examination, work of a similar nature with 
two-phase hydrogen (ref. 14) is used. 


SYMBOLS 


A area 

a thermal expansion coefficient 

B constant, 1.26 

C sonic velocity 

c specific heat 
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D droplet diameter 

d pipe diameter 

f frequency, Hz 

g gravitational constant 

h enthalpy 

hjg heat of vaporization 

J mechanical equivalent of heat 

K flow coefficient 

k thermal conductance 

L error producing variable 



p pressure difference 

Q energy added 

q energy loss 

R V r g 

r radius 

T temperature 

T q calibration temperature 

V velocity 

V free-stream velocity 
s 

v specific volume 

w weight-flow rate 

AX sampling error 

x quality 
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x free -stream quality 

Y expansion factor 

y moisture content, 1 - x 

a thermal diffusivity, k/pCp 

/3 orifice-pipe diameter ratio 

ju fluid viscosity 

v kinematic viscosity 

p density 

r time constant 

co frequency, rad/sec 

Subscripts: 

f saturated liquid or droplet 

g saturated vapor 

is isentropic 

o near measurement point or orifice 

x mixture 

1, 2, 3 measurement locations 


PRELIMINARY CONSIDERATIONS 
Description of Fluid 

In the quality range of interest, 0. 85 to 1. 0, the fluid can be described as a mist; on 

a volumetric basis, it is mostly vapor with the liquid phase dispersed in the form of 

small droplets . The droplet diameters are estimated to range between 1 and 100 microns 

2 8 

at concentrations between 10 and 10 droplets per cubic centimeter, depending on the 
method of vapor generation, phase densities, stream quality, and stream turbulence. 

For typical potassium and sodium systems, the Reynolds numbers range from 10^ up- 
ward, which indicates turbulent flow conditions. Because of the high droplet concentra- 
tion and turbulent flow conditions, the fluid is considered to be a continuous and homoge- 
neous mixture of the two phases. 


4 



Evaluation of Thermodynamic Properties 


if 


In the measurement of quality, thermodynamic properties must be used. The liter- 
ature reveals a number of sources for property data. A comparison of data (ref. 15) , 
however, indicates that significant differences exist among the property values given in 
various sources. As an example, an 8 -per cent variation in values for the heat of vapor- 
ization of potassium at 1800° F can be calculated from figure 13 of reference 15. These 
differences result in systematic uncertainties in the measurement. They will not enter 
into the random error calculations presented for the calorimetric and flow measurement 
techniques. In the section, Effects of property variation on quality computations, the 
difference in quality values resulting from the use of two tables of thermodynamic prop- 
erties is given. 

Thermodynamic properties for use in quality determinations are obtained from 
tables or curves with the use of the primary measurements of temperature and/or pres- 
sure as indices. The accuracy by which the properties can be found is a function of the 
accuracy involved in the primary measurement and may be determined by standard error 
analysis methods. To evaluate the uncertainty in a property requires knowledge of the 
partial derivative of the property with respect to temperature and to pressure and esti- 
mates for the errors in these measurements. The partial derivatives are computed on a 
finite difference basis from property tables. The choice of a particular table is not crit- 
ical since reference 15 shows these partial derivatives to be similar in the temperature 
range of interest. 

The errors in temperature and pressure measurements are considered random. 

The probable error of an absolute pressure measurement in liquid metals is estimated at 
±1 percent of full scale or ±0. 1 pound per square inch, whichever is larger. The as- 
sumed probable error in a temperature measurement is ±3/8 percent of the reading. 
These estimates are considered optimistic but are useful in establishing the potential 
value of the techniques considered. 

A summary of/the uncertainty in the saturated thermodynamic properties of potas- 
sium at 1600° and / 2200° F is given in table I. The computations were based on data 
from reference 16. The properties are shown in two columns depending on whether the 
temperature or the pressure was the measured variable. Note that the largest uncer- 
tainty involves the determination of the specific volume of the saturated vapor. 
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TABLE I. - ERROR IN PROPERTY 


EVALUATION FOR POTASSIUM 


[AP/P, 0.01; AT/T, 0.00375.] 


Error producing 
variable, 

L 

Saturation temperature, 

°F 

1600 

2200 j 



Error, 

percent 



i^AP 
l ap 

— — AT 
L 3T 

I^AP 
l ap 

I9L AT 
L 3T 

Specific volume 





Liquid, v f 

0.053 

0. 13 

0. 10 

0.21 

Vapor, v g 

.83 

2.2 

. 87 

1. 8 

Enthalpy 





Liquid, h f 

0. 11 

0. 30 

0. 17 

0. 35 

Vapor, h g 

.088 

.021 

.018 

.036 

Heat of vapori- 

0.045 

0. 12 

0.093 

0. 19 

zation, h y 






ANALYSIS OF MEASUREMENT TECHNIQUES 
Calorimetry 

Calorimetric quality measurement techniques utilize a steady-flow process in which 
the two-phase fluid is either heated, condensed, or expanded. The process continues 
until a single-phase fluid of measurable energy content is achieved. By measuring the 
amount of energy required to achieve the final state, the initial energy content and there- 
fore quality of the two-phase fluid is determined. 

A schematic diagram of a calorimetric system is shown in figure 1. The system 
consists of a sampling probe, a process chamber, and a flowmeter. The fluid enters the 
system at station 1, where a temperature or pressure measurement establishes the ini- 
tial saturated state. Energy is removed or added, or the fluid is expanded until a single - 
phase fluid is achieved at station 2. The weight flow rate of the fluid is then measured. 
The steady-flow energy equation for the process, with gravitational effects neglected, is 



2Jg 


. Q 
+ h. + = 

w 




( 1 ) 
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( 2 ) 


The enthalpy h 1 is related to quality by 

h, = h g> 1 - (1 - x)h tg _ j 

Combining equations (1) and (2) and neglecting the kinetic energy terms result in the fol- 
lowing expression for quality: 


x=l- V 1 '" 2 -— Q— t — 3— (3) 

hfg, 1 wh fg, 1 wh fg, 1 

For systems in which the velocities are in the range of 100 feet per second or less, the 
kinetic energy term is small in comparison with the other terms and can be neglected. 

If stream velocities are much higher or uncertain in magnitude, they must be considered 
in an error analysis. 

Equation (3) includes two energy terms. The first term Q/wh f g j represents the 
measured amount of energy required to achieve a single-phase fluid at station 2. In this 
derivation, the term is positive for heating, negative for condensing, and zero for thrott- 
ling calorimeters. The second term Q/wh^ ^ considers heat losses from the process 
chamber that are functions of system design and environmental operating conditions. 

The actual value for the heat loss may be determined experimentally and a systematic 
correction applied to the measurement. This term is treated as a constant in the error 
equations. 

The root -sum -square error as discussed in this report is defined by the following 
expression: 


Ax 


2 




( 4 ) 


Simplified expressions for the root-sum-square error in each of the calorimetric 
systems can be obtained by making use of the following assumptions and approximations: 
(1) Negligible superheating or subcooling occurs in the process chamber. This as- 
sumption implies that (h^ ^ - hg) /h^ ^ is 0. 0 for a heating calorimeter and is 1. 0 for 
a condensing calorimeter. 
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(2) The following approximations are made: 


dh g, 1 _ dh 2 _ ^fg, 1 
h g, 1 h 2 h fg, 1 

(3) The heat loss correction q/wh fg ^ is assumed to be 0. 01 with error of ±10 per- 
cent, and is assumed to be constant. 

Using the values for potassium at 2200° F from reference 17 produces the following 
error equations: 

Heating calorimeter: 



(5a) 


Condensing calorimeter: 



(5b) 


Throttling calorimeter: 
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In the following paragraphs, various physical configurations of the calorimeters are 
described. Based on these configurations, estimates for the error terms in equations (5) 
are made. Also described are some of the operational characteristics peculiar to each 
calorimeter. 

Heating and condensing calorimeters . - For the heating calorimeter, the most con- 
venient source of energy input to the process chamber is an electric heater. The prob- 
able error in measuring the heat input is 1/2 percent. A condensing calorimeter re- 
quires a heat exchanger in which the actual energy measurement is probably made on a 
second fluid. The probable error of such an energy measurement is estimated at 1 per- 
cent. Most calorimeters operate with flow rates in the range of 20 pounds per hour. 

For this range, thermal flowmeters are practical, with a root-sum-square error of 
1^ percent. Some design and operating problems for thermal flowmeters include mini- 
mizing heat losses and measuring the fluid bulk temperature. 

For a condensing calorimeter, a volumetric flowmeter might also be considered. 
This device involves the measurement of a time interval required to fill a known volume 
with liquid. The uncertainty in such a measurement system would be comparable to that 
achieved in the thermal devices. The magnetic flowmeter does not appear applicable to 
this problem because of the small flow rates involved. 

The principle of a calorimeter involves achieving a single-phase fluid at the outlet. 
Thus, ideally, the energy content is directly measurable. In a heating calorimeter, 
however, liquid carryover may occur; thus, some uncertainty is introduced in the state 
of the fluid. The extent of this problem depends on initial droplet size, fluid turbulence, 
surface-tension effects, and the amount of droplet coalescence. In contrast, the con- 
densing calorimeter should be capable of operating such that only liquid is present at the 
outlet, and therefore, the question of phase composition is not involved. 

Throttling calorimeter. - The throttling calorimeter utilizes a steady-flow adiabatic 
expansion process in which no work is done on or by the fluid. Consequently, the only 
error involved in this device is that caused by determining enthalpies. A discussion of 
throttling-calorimeter feasibility is given in reference 12. The report indicates that 
the throttling calorimeter is limited from considerations of thermodynamics to measure- 
ments in the 0. 91 to 1. 00 quality range for potassium. The limitation is based on inlet 
temperatures between 1600° and 2200° F and on expansion pressures of 1 to 10 pounds 
per square inch absolute in the process chamber. 

From equation (5c), the uncertainty in a measurement of quality approaches 0. 0025 
with a throttling calorimeter based on Ahjg l/^fg l = 001. limit of error was es- 

timated at 0. 01 in ref. 12, as opposed to the root-sum-square error used in this re- 
port. ) Because of the simplicity of a throttling calorimeter and from a more thorough 
heat-transfer analysis, the heat loss error was estimated at 0. 005 in reference 12 and 
was not included as a correction in the measurement. 



Summary of errors in calorimeters . - A summary of the root -sum -square errors 
for calorimeters is given in table n. Indicated are the assumed values of error for the 
latent heat, weight-flow, and energy measurements. Under these conditions, less error 
is involved in using throttling and heating calorimeters than in using condensing calo- 
rimeters. For covering a wide range of qualities, the use of a heating calorimeter 
should be considered. 


TABLE H. - ROOT-SUM-SQUARE ERROR IN CALORIMETERS 


Calorimeter 

Error term 

Quality 

Ah f g 

Aw 

AQ 

Q 

q 

Aq/whfg 

0.85 

0.95 

V 

w 

wh f g 

q/wh fg 

Root - sum - square 
error 

Heating 

0.001 

0.0125 

0. 005 

0.01 

0. 10 

0. 0034 

0. 0028 

Condensing 

.001 

.0125 

.01 

.01 

. 10 

.014 

.015 

Throttling 

.001 



.01 

.10 


. 0027 


Sampling error . - Obtaining a representative sample of a two-phase fluid is a prob- 
lem and is an important prerequisite for achieving reliable calorimetric quality mea- 
surements. Even when the liquid and vapor phases are homogeneously dispersed in the 
flow line, as assumed herein, an error caused by sampling can arise. This error re- 
sults from the difference in momentum of the liquid drop and the vapor. A number of 
reports are available (refs. 7, 9, and 12) that recommend methods for minimizing sam- 
pling error. One basic recommendation is that the probe and flow system be designed to 
ensure that the velocity at the probe inlet is equal to that of the free stream at the mea- 
surement plane. In reference 12, the effects of velocity differences on sampling error 
are discussed, and the following equation is derived, which relates sampling error to the 
velocity ratio and quality: 


Ax = x„ 
s 



where 


V = 
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For an allowable sampling error, the required degree of control and knowledge of the ve- 
locities can be obtained by solving the equation for V. 


I 



( 6 ) 


As an example, if the allowable sampling error were 0. 01, the velocity at the inlet 
must be controlled to within 20 percent of V_ when measuring qualities in the 0. 95 
range. However, for measurements in the 0. 85 range, the control must be to within 
8 percent. Thus, to minimize the sampling error in the lower quality range, precise 
control of free-stream and sampling-probe velocities is required. In many liquid-metal 
applications it is difficult to control the velocities to the required degree, and, there- 
fore, an analysis of other measurement techniques that do not necessarily require sam- 
pling was undertaken. 


Flow Measurement 

In reference 18, Murdock derives an equation relating the total weight flow of a two- 
phase fluid passing through a sharp-edged orifice to the measured pressure drop, the 
moisture content, and fluid metering constants. This equation was developed from a ra- 
tional approach to two-phase flow measurements and contains an empirical constant de- 
termined from correlating many two-phase flow measurements. The equation was de- 
veloped for use in steam -water flows in large power plants and for metering in natural 
gas piping systems. Since the moisture content is equal to 1 - x, this metering system, 
in conjunction with a second measurement of total weight flow, could determine quality. 
The potential advantages of this method are that measurement of the entire fluid flowing 
in a pipe eliminates sampling and that absence of heat losses eliminates a source of cor- 
rection required for calorimetric devices. Thus, the evaluation in some detail of the 
applicability of this flow-measurement technique for liquid-metal quality determination 
is of interest. 

The total weight-flow rate of a two-phase fluid passing through a sharp-edged orifice 
is found from the following equation developed in reference 18: 
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w = 


CO 



The constants in equation (8) include flow coefficients, the orifice area, and an ex- 
pansion factor. In fluid metering, the flow coefficients are determined with the use of a 
Reynolds number correlation. For liquid-metal applications, both air and water may be 
used to obtain the necessary Reynolds number range for determining the constants Kg. 
and Kj. The fluid expansion factor takes into account the compressibility of the fluid 
and is obtained from tables by using the ratios of specific heats and measured pressures 
as indices. 

Error evaluation for two-phase flowmeter technique . - For the purposes of this sec- 
tion, all sources of error in the measurement of quality with a flowmeter are considered 
as random errors. The total error in the measurement is then taken as the root sum 
square of the individual errors (eq. (4)) , which can be written as 
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where 


K Y A / \l/2 

M - K g Y g A Y2Gp\ 

W VW 


k /v f y / 2 

N = B — ^ Y (— ] 

K f \g/ 


The error terms are evaluated in appendix A. A summary of the terms is as fol- 


lows: 


AK 


K 


2. lXlCf 3 


g 


AY 


I = 4xl0 -3 


g 


I /_l\ = 5xl0 -3 


2 \ v 


g. 


2.5X10" 3 


= 10xl0~ 3 
w 
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A /k P\ _ 5xlO -3 
2 \P / 


AB 

B 


5X10" 3 



Basing the values of M and N on the specific volumes for saturated potassium at 
1600° F (ref. 16) and assuming that = 0. 97, = K^, and x = 0. 90 produces repre- 

sentative values for M and N of 0. 90 and 0. 052, respectively. Thus, the root-sum- 
square error as given by equation (9) is 0. 014 for the flow measurement technique. This 
error would be equivalent to that computed for either the throttling or heating calorim- 
eter when sampling errors (0. 01) are considered in these measurements. 

Differential pressure measurement. - The differential pressure measurement is the 
primary measurement required when the flow measurement technique is used. A gen- 
eral indication of the differential pressure range can be obtained by considering the 
pressure drop developed in a sharp-edged orifice through which a saturated liquid-metal 
vapor is flowing. Figure 2 shows the pressure difference p as a function of fluid veloc- 
ity with the orifice coefficient assumed to be 0. 6. Lines are shown for potassium at tem- 
peratures of 1600° and 2200° F. For practical meter systems, the throat velocities 
would be about 500 feet per second, which corresponds to differential pressures between 
6 and 30 pounds per square inch. For a two-phase fluid, the differential pressure is 
somewhat higher because it is approximately inversely proportional to quality for a given 
fluid velocity. 

The measurement of low pressures in liquid-metal systems is difficult. Fortu- 
nately, there are a few methods, such as inert gas bleed systems, which are capable of 
achieving differential pressure measurements (refs. 19 and 20) accurate to within 
±0. 1 pound per square inch. Also, commercial high-temperature differential pressure 
transducers are available, and if care is used in their installation and operation, accu- 
racies approaching 1 percent of the reading (as estimated in the error analysis) are ob- 
tainable. Flowmeter systems such as these have one limitation; namely, a given system 
is capable of making accurate measurements only over a small range of weight flow 
rates around the design point. This limitation occurs because the weight flow rate is 
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proportional to the square root of the pressure difference p. Consequently, for a given 
uncertainty in p, the error becomes large at flow rates less than the design point. Ex- 
amination of equation (8) and the coefficients M and N leads to a similar conclusion for 
the measurement of quality. This particular characteristic may be a deciding factor in 
the choice of the flowmeter technique for quality measurements. 

Effects of property variation on quality computations . - A comparison of the effects 
of differences in thermodynamic property data on the computed values of quality for the 
calorimetric and flow measurement techniques is now made. In appendix B are derived 
equations (Bl), (B2), (B3), and (B5), which relate the difference between the values of 
quality computed from the data of reference 16 and 21 to inlet saturation temperature. 
The derivations were based on the assumption of no superheat or subcooling at the outlet. 
Also, use was made of linearized relations between the saturation properties and tem- 
perature. From the equations, the following conclusions are reached: For the two sets 
of property data used, the difference in computed values of quality for the heating calo- 
rimeter approaches 6 percent of the moisture present in the stream while that for the 
condensing calorimeter approaches 6 percent of the measured value of quality. For the 
throttling calorimeter, the differences approach one-half the amount of moisture present 
in the stream while the difference for the flowmeter system ranges between ±3 percent of 
the measured value of quality. 


Sound Velocity 

The sonic velocity in a two-phase fluid is theoretically a function of fluid quality. 

The measurement of sound velocity, therefore, might be used as the basis for a quality- 
measuring system. Such a system would consist of a sound generator and an acoustic 
transducer located a known distance from the generator. The advantages of this system 
include freedom from sampling errors and heat loss corrections, provided that the re- 
quired instrumentation can be located in the stream. In the following sections, an analy- 
sis of the thermodynamic process is given along with a discussion of the heat transfer 
and mechanics of the vapor and droplet interaction. 

Process equation . - The thermodynamic equation relating sonic velocity to fluid 
properties for any fluid is 


C 


2 



( 10 ) 
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TABLE IB. - SONIC VELOCITY 
OF POTASSIUM 


where 


Temperature, 

°F 


Quality, 

X 


0.85 

0.90 

0.95 

0.98 

0.99 

Sonic velocity, ft/se< 


1600 

1558 

1606 

1653 

1681 

1690 

1700 

1586 

1636 

1684 

1713 

1722 

1800 

1612 

1663 

1713 

1742 

1751 

1900 

1635 

1688 

1739 

1769 

1778 

2000 

1657 

1710 

1763 

1793 

1803 


H = 

W/.„ 


Papers have been published giving the sonic 
velocity as a function of quality and tempera- 
ture (e. g. , refs. 1 and 14) for many fluids. 
Table HI gives values for the sonic velocity 
of potassium that were computed on a finite 
difference basis with the use of equation (10) 
and the equations of state as given in refer- 
ence 17. The table indicates that the sonic velocity changes approximately 10 feet per 
second for a 0. 01 change in quality. If the sonic velocity can be measured to within 
1/2 percent, the resulting error is approximately 0. 008 in quality. Also, this table 
shows that the effects of temperature on the sonic velocity are not great, with a change 
of 30 feet per second for a 100° F change in temperature. With a probable error of ±3/8 
percent for temperature measurements, an additional uncertainty of 0. 0018 (at 1600° F) 
is introduced into the measurement. The resulting root- sum -square error is 0. 008. 

When use is made of equation (10), the fluid must be in thermodynamic equilibrium 
during the passage of the acoustic wave. To maintain equilibrium, both phases must be 
at the same temperature and move as a unit at all times. Since it is obvious that these 
conditions cannot be achieved over wide ranges of acoustic frequencies, it is of interest 
to examine the mechanical and thermal aspects of this problem. 

Heat-transfer considerations. - Consider a mixture in which both phases of the fluid 
are at the same saturated conditions. On a typical temperature-entropy diagram (fig. 3) 
the mixture is assumed to exist at state point 1. The propagation of an acoustic pulse 
will result in a compression and a rarefaction of the fluid. For the isentropic compres- 
sion process, state point 2 may be achieved. This process results in a change in the 
quality of the mixture and, therefore, in mass and heat transfer between phases. 

Since energy cannot be transferred instantaneously, a relaxation time is involved in 
the process. For the purposes of establishing the relaxation time, the mixture is as- 
sumed to consist of a number of spheres equal to the number of droplets present. Each 
sphere consists of a liquid droplet surrounded by a vapor shell. The volume of the shell 
is so proportioned that the quality of a given sphere is the same as that of the mixture. 

To simplify the problem further, it is assumed that no heat is transferred between the 
various spheres. The relaxation times or time constants for the process of interest are 
those which are characteristic of the temperature distribution in a droplet or shell when 
a step change in the interface temperature occurs. The equations for these time con- 
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TABLE IV. - RELAXATION TIME CONSTANTS 


FOR POTASSIUM AT 2200° F 


Quality, 


Droplet diameter, 

fim 

X 

2 

10 

50 

100 


Time constant of vapor shells, r g , sec 

0.80 

2xlO" 5 

5X10' 4 

lxlO' 2 

5X10' 2 

.90 

4x10”® 

lxlO” 3 

3xl0” 2 

lxlO" 1 

.95 

1X10” 4 

2xl0” 3 

6xl0’ 2 

2X10” 1 

.99 

5xl0” 4 

lxlO" 2 

3x10” 1 

1. 2x10° 


Time constant of liquid drops, r, sec 



2xlO -9 

6X10" 8 

2xl0” 6 

6X10” 6 


stants are taken from reference 14. For a 
sphere, the relaxation time constant is 
found from 


T f = 


(H) 




The relaxation time constant for a spheri- 
cal shell is 


T g = 


l f 


(r - 2 - l) 


( 12 ) 


a 6 
g 


where R is the ratio of the droplet radius 
phase density by 


R = 


to sphere radius and is related to quality and 

1/3 


1 + 


x \/Pf\ 

1 - X/ Miff) 


(Eq. 3. 5 in ref. 14 differs somewhat from this equation. ) The angle 6 is found from the 
following transcendental equation: 


tan 6 

Relaxation time constants for potassium at 2200° F are given in table IV. The table 
indicates that the time constants for the size of droplets expected in liquid- metal systems 
are very short, less than 10” ^ second. On the other hand, the time constants for the va- 
por shells approach seconds. Because of these long time constants, slow pressurization 
rates and, therefore, low acoustic frequencies will be required to maintain equilibrium. 
Note that the foregoing calculations do not include the heat required to effect a change in 
quality; thus, the actual time constants would be longer. 

Mechanical considerations . - A second limitation on the operating frequency of an 
acoustic system results from the requirement that the vapor and droplets move as a 
unit. In reference 14, general equations are given for the passage of an acoustic wave in 
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a fluid containing particles and vapor. These equations were derived by considering the 
viscous forces acting on the particles and their resulting motion as a function of the va- 
por motion. It was assumed that the vapor and droplets move in a sinusoidal manner and 
that the movement of the droplet lags that of the vapor by some phase angle e. The fol- 
lowing equation from reference 14 gives the ratio of the droplet velocity to that of the va- 
por. The droplets were assumed to be rigid spheres. 



(13) 


where 


$ = 




P 


g 


P f 


4 /3r f y pr f J p { 


(ft 


$ + 1 

Values for the velocity ratio are shown in figure 4 as a function of a diameter- 
frequency parameter fD^ for potassium at 1600° and 2200° F, where (/3r^)^ = 

These curves show that the droplets will move with the vapor for values of fD^ < 10^. 
Thus, for a 50 -micron droplet, the system must operate at 4 hertz or less to minimize 
the relative velocity. This operating frequency also approaches that necessary for 
maintaining temperature equilibrium as previously discussed. 

The effect of relative velocity between the vapor and liquid droplets on the sonic ve- 
locity of the fluid is shown in figure 5. The curves were computed from an equation 
which is a solution to the wave equation for a pressure pulse transmitted through a vapor 
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in which are dispersed many particles. The equation is derived in reference 14 and is 
as follows: 




4>(1 + <l?) + 0 2 


P s 

l - -E 


V x ) 

V 

_(i + $) 2 + e 2 _ 


p f. 



11 - 1/2 


(14) 


2 4 

For fD < 10 , the sound velocity is constant and a function of the mixture quality. As 
fD 2 increases, the presence of the droplets affects the sound velocity to a lesser degree. 

This examination of the heat transfer and mechanics of the phase interaction indicate 
that the use of a sound velocity measurement for quality determination requires operating 

4 

at a value of the diameter-frequency parameter less than 10 . In the section, Descrip- 
tion of Fluid, the droplet diameters were estimated to be between 1 and 100 microns. 
Covering this range of diameters and maintaining a diameter-frequency parameter value 
of 10 4 would require measuring frequencies approaching 1 hertz. In practice, liquid- 
metal systems often exhibit pressure fluctuations at these low frequencies. Therefore, 
signal to noise ratio problems would be encountered, which introduce uncertainty in a 
measurement. Since quality measurement techniques cannot be calibrated to overcome 
such problems, the sonic velocity measurement technique does not appear practical. 
However, in determining the presence of small droplets, the technique may be of some 
use since higher operating frequencies may be utilized and the signal to noise ratio 
thereby maximized. 


CONCLUSIONS 

In an analysis of three liquid-metal quality measurement techniques, the use of cal- 
orimetry was concluded to be the most practical method available at this time. On the 
basis of applying a systematic heat loss correction, estimated root-sum-square errors 
for these devices were 0. 015 for a condensing calorimeter, 0. 0028 for heating calorim- 
eters, and 0. 0027 for throttling calorimeters. The major source of error in these de- 
vices resulted from errors in the weight flow measurement. Also, sampling of the two- 
phase stream presents a problem when using calorimetric methods with sampling errors 
approaching 0. 01. 

Two-phase flow measurement techniques have a potential root- sum- square error of 
0. 014. The measurement of the small differential pressures is a problem in the use of 
this technique. The technique is limited to measurements near design conditions be- 
cause of the functional relation between quality and differential pressure. 
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The use of a sonic velocity measurement technique does not appear practical at this 
time. The main problem is associated with achieving thermodynamic equilibrium during 
the passage of a pressure wave. Because of the required diameter-frequency parameter 
range for achieving equilibrium, the environmental conditions, and uncertainty in droplet 
sizes, this technique is limited in application. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 20, 1967, 

120-27-02-04-22. 
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APPENDIX A 

FLOWMETER ERROR TERM EVALUATION 
Flow Coefficient 

The flow coefficients AKg/Kg and AKj/'K^ are a function of Reynolds number and 
are usually obtained from a calibration with air or water used as the working fluids. The 
accuracy of a calibration is in the 1/8 percent range. The total error in knowledge of 
the flow coefficient is found in the following manner: Assume that the functional relation 
exists such that 


K = K + f(Re) 

g g, o v ' 


where Re is the Reynolds number. Let Kg Q » f(Re). Then 


AK 


K 


g_ 


r AK \2 

-1^2 \ + _J ARe 


K 


g,o 


K_ „ 3Re 
g, o 


1/2 


The first term in equation (Al) is taken as the calibration error: 


AK 


0.00125 


K 


g,o 


(Al) 


The partial derivative in equation (Al) will be evaluated on a finite difference basis from 
tables in reference 22. The sharp-edged orifice under consideration will be for a 
2-inch pipe with flange taps, Reynolds number of 10 000, and an orifice-pipe diameter 
ratio 0 of 0. 25. For this configuration K = 0. 61, and from table 6 of reference 22, 

_7 

the partial derivative 3f/3Re is found as 9X10 . With Re = pVD/p, the error in the 

Reynolds number is 


ARe 

Re 
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The error in each of these terms is assumed as follows: 


*£= 0.01 

P 

AV 


V 

AD 

D 


= 0. 05 


= 0. 005 


0. 10 
M 

and the error in Reynolds number becomes 

0.112 

Re 

where the major portion of the error is involved in the viscosity. Substitution of the 
specified values for the terms in equation (Al) leads to a total error in the flow coeffi- 
cient AK /K of 2. 1X10' 3 . 
g / g 


Expansion Factor 

The expansion factor AY g /Y g for a fluid passing through an orifice is a function of 
the pressure ratio and adiabatic expansion coefficient: 

Y = f(r’,y) 


where 

r’ pressure ratio, 
y adiabatic expansion coefficient 
Pj upstream pressure 

The root-sum-square error in the expansion factor is 
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VV k 


(A2) 


The partial derivatives can tjre evaluated from figure 72 of reference 22 on a finite differ- 
ence basis. For values of k = 1. 4, r* = 0. 9, Y = 0. 97, and 0 = 0. 25, the partial deriv- 
atives are 


0r’ 

M. = o. 03 
0k 

If the pressures are assumed to be measured to 1 percent and Ak = 0. 02, the total error 
AY /y in the expansion factor computed from equation (A2) is 4xi0 -3 . 

o & 


Orifice Area 

The total error A A/A in knowledge of the orifice area must take into account the 
area change with temperature. The area of the orifice is given by 

A = 4 K (1 + aAT) ] 2 


where 


d„ diameter measured at T 
o o 

The total error in the area determination is given by 


AA 


= 2 < 


r 

J 

\h\\ 

' a(T - T q ) 

2 

/Aa\ 2 , / a < T - T o>\ 2 


v_ 

YV 

1 + a(T - T q ) 


U / V T - T o) 

r 


1/2 


If the values for the terms, Ad o/ /d Q = 0. 001, T - T q = 1500, a = 10 -5 , Aa/a = 0. 05, and 
A(T - T q )/T - T q = 0. 00375 are assumed, the total error AA/A in the area is 

2. 5X10“ 3 . 
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Factor B 


The uncertainty in measuring total two-phase weight flow with the use of equation (7) 
was given as ±1. 5 percent in reference 18. A portion of this is involved with uncertainty 
in the empirical correlating factor B. Reference 18 indicates that the standard deviation 
of B was 0. 75 percent. Therefore the random error in the factor B is assumed equal 
to its probable error (which is two-thirds of the standard deviation) or AB/B = 0. 005. 

The uncertainty in the remaining factors in equation (9) are assumed as follows: the 
weight flow Aw/w was 0. 01; the differential pressure measurement Ap/p was 0.01; 
and the specific volume measurements Av^/v^ and Av^/v^ were 0. 01. Thus, values 
for all terms within the summation signs are established. 
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APPENDIX B 


DERIVATION OF DIFFERENCES 
Calorimeters 

Under the assumption that no subcooling or superheating occurs in the devices, the 
enthalpy hg in equation (3) becomes either a saturated vapor enthalpy or a saturated liq- 
uid enthalpy for calorimetric devices. The differences in the values of quality computed 
from two sets of property data then follow directly from equation (3) as 

Heating calorimeter: 

ix HC = 

Condensing calorimeter: 

Ax CC 


h fg( h g,l~ h g,2)- h fg( h g,l~ h g,2) 

bfg h fg 

For this appendix, the primed variables represent data from reference 16, while the un- 
primed variables are data for potassium from reference 21. To eliminate the enthalpy 
terms as variables and obtain expressions for these differences as a function of satura- 
tion temperature, plots of hj g against h fg and h^ against h g were constructed. 

The points were plotted for specific values of saturation temperatures, and from the 
curves the following linear relations were obtained: 

h, = 0.773 ht + 212 
fg fg 

h = 1. 75 h* - 882 
g g 

Also, from a plot of hj g against temperature, the following relation was obtained: 


Throttling calorimeter: 

Ax TC = 
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h* g = -0. 152 T + 1039 


Eliminating the enthalpies in the difference equation and using equation (3) yielded the 
following difference equations: 


Ax Hr =(l-x»)f-gi-- 0 - 0345 T 
\1015 - 0. 117 T; 


Ax rr = -x* (A 4 - 0-0345 T> 
^ \1015 - 0. 117 Ty 


Ax, 


TC 


= (1 - x ') (- 806 ~ °- 149 T 
\1015 - 0. 117 T 


(Bl) 

(B2) 

(B3) 


Flow Measurement System 

Variations in the available values of the specific volumes of the saturated fluids re- 
sult in differences in computed values of quality in the flow measurement technique. To 
obtain the difference equation, the following substitutions are made for terms in equa- 
tion (8): 


C 3 = ~y Y y A (2Gp) 1/2 

w 


K„ 

C A = B -i Y„ 
4 K, S 


Z = 



u = 



The terms of interest are Z and U. The equation for the flow measurement technique 
becomes 
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Let the primed variables Z' and U' refer to specific volumes given by reference 16 
and the unprimed variables Z and U refer to the data for potassium from reference 21. 
The difference equation becomes 


Ax 




1 - C 4 U 


where Cg has been eliminated by use of equation (B4). 

The terms Z/Z’, U’ Z/Z *, and U were plotted as a function of saturation tempera- 
ture and the following relations were obtained: 


— = 1. 177 - 9. 11X10" 5 T 
Z’ 


U’ — = 9. 56X10” ^ T - 0. 108 
Z* 


U = 9. 56X10” 5 T - 0. 109 


The values for the terms in the constant C 4 were assumed as 

B = 1. 26 

K f= K g 

Y g =1 -° 

C 4 = 1. 26 

The difference equation for the flow- measurement technique becomes 


Ay _ v ,/9. 09X10”° T - 0. 175 

ax FM “ x \ 

a. 13 - 12.0X10"° T 


(B5) 


where 1. 31 x* » 1. 26xl0" 3 . 
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Throat velocity, ft/ sec 



Figure 1. - Calorimeter. Energy added Q is positive for heating 
calorimeter, negative for condensing calorimeter, and zero 
for throttling calorimeter. 
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Figure 2. - Variation of differential pressure with velocity. 




Figure 3. - Temperature -entropy diagram. 



Figure 4. - Variation of velocity amplitude ratio with diameter frequency 
parameter. 
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Figure 5. - Variation of sound velocity ratio with diameter -frequency 
parameter for potassium at 1600° F. 
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